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Цель работы состоит в изучении корректирующей способности 
нового мягкого декодера кодов Рида — Маллера. Метод до- 
стижения цели заключается в экспериментальном исследова- 





нии декодера с использованием специально построенной ИМИ- 
тационной модели помехоустойчивого канала передачи дан- 


ных. Источник и приемник сообщений модели оперируют 


цифровыми данными, заданными полем ЁРз. Линия связи по- 
строенной модели в зависимости от настроек выдает цифровые 
или непрерывные сигналы. В случае непрерывных сигналов 
рассматриваются два варианта базовых искажений сигнала и 
их комбинации. Помехоустойчивость моделируемых каналов 
связи обеспечивается применением кодов Рида — Маллера 
второго порядка, заданных над полем Ез, и нового декодера 
мягких решений для этих кодов. Результаты проведенных ими- 
тационных экспериментов показали, что исследуемый декодер 
как в цифровых, так и в полунепрерывных каналах позволяет 
исправлять ошибок больше, чем гарантируется минимальным 
кодовым расстоянием. Наибольтую эффективность декодер 
показал при использовании его в полунепрерывных каналах 
связи. Корректирующая способность декодера значительно 
зависит от типа линии связи и вида искажений, поражающих 
сигналы, и не чувствительна к местоположению ошибок внут- 
ри кодового слова. Сделаны выводы о возможности использо- 
вания нового декодера в каналах связи низкого качества для 
обеспечения помехоустойчивости, а также в ряде криптогра- 
фических приложений. 


Ключевые слова: троичный канал связи, троичные коды Ри- 
да — Маллера, декодер мягких решений, математическая мо- 
дель канала связи, экспериментальное исследование корректи- 
рующей способности кода, исправление ошибок за границей 
половины минимального расстояния кода. 


Введение. В разнообразных системах передачи и хранения информации для ее защиты от искажений ИСПОЛЬЗУЮТСЯ 


алгебраические помехоустойчивые КОДЫ [1-3]. Одной из акгуальных задач в Настоящее время является разработка ДЛЯ 
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= ннианый 
известных кодов новых декодеров, которые обладают какими-либо преимуществами по сравнению с известными де- 
кодерами. К таким преимуществам может относиться, например, увеличение скорости работы декодера или повыше- 


НИЕ корректирующей способности кодека. Так, в последние Годы Для обеспечения помехоустойчивости актуально 





применение декодеров мягких решений (ДМР), особенность которых состоит в том, что принятые из канала данные 
вводятся в декодер, минуя демодулятор [1]. Обычно использование ДМР дает лучшие результаты по сравнению с де- 
кодированием жестких решений, когда на вход декодера поступают значения с выхода демодулятора, преобразующе- 
го данные из канала в слова над фиксированным конечным алфавитом. Эффективность ДМР основана на том, что в 
отсутствие демодулятора не накапливаются оптибки квантования, однако обычно декодеры с технологией ДМР обла- 
дают большей сложностью — например, [2]. 

В работе [4] построен мягкий декодер кодов Рида — Маллера второго порядка, заданных над полем Галуа Ё3з 
(далее АМ3(2, т), т — параметр кода). При этом за основу взят декодер двоичных кодов Рида — Маллера второго 
порядка В. М. Сидельникова и А. С. Першакова |5], обладающего значительной корректирующей способностью (он 
исследовался в работах [6—7|). Однако нет теоретической или экспериментальной оценки корректирующей способно- 
сти нового алгоритма. 
Цель данной работы состоит в экспериментальном исследовании корректирующей способности нового мягкого де- 
кодера [4] троичных кодов Рида — Маллера второго порядка при различных условиях его эксплуатации. Для дости- 
жения цели необходимо решить две задачи. Во-первых, построить модель канала с троичным входом. Источник и 
приемник сообщений данной модели оперируют цифровыми данными, заданными над полем Ё3з, а линия связи в зави- 
симости от настроек выдает цифровые или непрерывные сигналы над полем комплексных чисел. Во-вторых, провести 
экспериментальное исследование корректирующей способности нового декодера при различных видах искажений 
кодовых слов. С этой целью используется программная реализация построенной модели троичного канала связи. 
Модель троичного канала передачи данных © использованием канального ДМР-кода Рида — Маллера второго 
порядка. Рассмотрим элементы модели канала с троичным входом и схему прохождения данных по модели (рис. 1). 


Источник сообщений выдает информационные векторы 
сы А 
т = (т, то....,т,}е РЗ, 
К а | | й Н | 2." 
где Ру — линейное К-мерное пространство, заданное над полем Галуа Ру. 
Затем в кодере канала эти векторы обрабатываются с использованием линейного блочного кода ЮМ3з(2, 17 ) 


Рида =—= Маллера второго порядка длины п и размерности К(< п) ‚ заданного над полем Е ь 
В —— И у хл 
Сформированные кодовые векторы Е = В поступают в передатчик, который служит интерфейсом К ЛИНИИ 
связи и преобразует векторы С Е №; в 


= (21, 2..5 би {= СЗ) 2 


где С — Цоле комплексных чисел, сигналы 2$, в ль принадлежат мультипликативной группе 


СР ху ху 
(в — *@ ы корнеи третьеи степени из единицы. 
4=0,1,2 
Преобразование аддитивной группы поля Аз в мульгипликативную группу Сз происходит с помощью есте- 


ственного изоморфизма ф: Е — С. ‚ который определяется по формуле 


ф(Л)=е 


Сформированные векторы 2 =(21. 22....2„) передатчик на физическом уровне отправляет в линию связи. 
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Физический аналог сигнала 2, можно получить, например, с помощью модуляции с непрерывной фазой [2]. Диа- 


грамма пространства таких сигналов иллюстрируется рис. 2. 
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Рис. |. Схема прохождения данных в моделируемом канале 





Рис 2. Диаграмма пространства сигналов 
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В Силу искажений. действующих в линии связи, на выходе из нее формируются СИМВОЛЫ ИЗ ПОЛЯ С. Будем 
р ыы = я 
рассматривать два оазовых вида искажении элементов вектора а. @. В ЛИНИИ СВЯЗИ — а именно, искажения По фазе 


И ПО амплитуде. Под искажением По фазе бу дем Понимать ф азовый СДВИГ сигнала 


2} по единичной окружности. Искажением по амплитуде будет смещение сигнала с единичной окружности По радиу- 


СУ. Предположим, ЧТо Под воздействием шума координаты вектора 2 подвергаются различным комбинациям базовых 


И — И \ ` 
искажении и формируется вектор 7 ЕС’. Графическая иллюстрация искажении сигнала представлена в левои верх- 


ней четверти рис. 5: Так, сигнал ег? показан черной Точкой, голубая точка соответствует его искажению По амплиту- 


де, зеленая — По фазе, а желтая — комбинации двух видов искажения. 


Искажение сигнала 7 Фильтрация сигналов, 

т 
по амплитуде. по фазе исказившихся за пределы 
И комбинация искажений ДОПУ стимой области 





Искажение при жестких 
решениях приемника 


Рис. 3. Примеры искажения сигнала и его фильтрации 


— и 
Вектор 2’ = (© поступает на вход приемника, который в зависимости от настроек может выдавать мягкие 


или жесткие решения о принягом сигнале. Так, в случае жестких решений приемник преобразует ВХОДНОЙ вектор 
— и — -7Я К 
=’ = С в вектор У — (У ры У, ) = (3 ‚ Независимо от типа искажении, поразивших вектор в ЛИНИИ СВЯЗИ. В этом 


случае может быть использован, Например, принцип решающих областей. Таким образом, в результате воздействия 


21, 27 
НЕ 


— Я т хл 127р т 
ОШИОКИ СИМВОЛ 25, =@е т С ав 7: может переити в один из двух других возможных символов: у, =е Ее С. ь 


Ь = [3  а= Ь. В построенной реализации модели такой переход осуществляется равновероятно. На рис. 3 показано, 


‚Чл. 
3 


что сигнал г - обозначенный черной ТОЧКОЙ, в результате искажения перешел в сигнал, обозначенный коричневой 
точкой. а приемник с жесткими решениями преобразовал его в сигнал еб р 

В случае мягких решений приемник фильтру ст амплитуду каждого сигнала 2 ‹ таким образом, чтобы значе- 
ния сигнала у, на выходе приемника принадлежали допустимой области 


веса 51/5}, 


А 


где 5(Е (0:1]) — параметр приемника. 
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= (у ‚) ==” Лействие фу т - 
У=(У.. ....У)ЕЯ.. Действие фильтра иллю 
стрируется на правой верхней четверти рис. 3: сигналы, обозначенные красными точками, попавшие в результате ис- 
кажений в линии связи за пределы допустимой области, переводятся приемником в сигналы, лежащие в области допу- 


стимых значений, с помощью смещения их по радиусу к ближайшей точке на границе области =.. Допустимая об- 
ласть выделена на рис. 3 серым цветом. Используя принятую в теории связи терминологию, можно говорить, что в 
случае работы приемника в режиме жестких решений реализуется цифровой канал, а в случае работы в режиме мяг- 
ких решений — полунепрерывный канал. 


Я 


Вектор у из пространства = ИЛИ С Гы =”) с выхода приемника В зависимости от режима заботы прием- 
3 5. 


с 
с 


г о —= А 
ника) направляется В декодер МЯГКИХ решении, цель которого —_ Восстановить информационный вектор те 5 ‚ ПО- 


сланный ранее источником сооощении. Результат декодирования у (Е 7: } поступает получателю сооощения. В зави- 


симости от уровня повреждения вектора 2 в канале связи результат декодирования может совпадать с исходным век- 
тором или отличаться от него. Если Ш =У ‚ то принято говорить о верном декодировании, иначе говорят об ошибке 
декодирования. 
Определение кодов Рида — Маллера ЮМз(2,т). Используя [4|, [8], определим троичный код Рида — Маллера 
ЮМз(2,т) второго порядка с параметром Ш> 2, следующим образом: 
— = 2 Я 
ВМ, (",т) = {( (о)... Л(а,))1 Л © ЕЗ[х,,...х„1 фе Е, 
эт 

где п=3 — длина кода, множество 


бы 0, } (5. - (0; о? м | 


представляет собой упорядочение всех точек векторного пространства Еу. 

Далее в работе используем следующее упорядочение: 

— по целочисленной сумме координат р(@) вектора ае Е" как натуральных чисел от меньшего к больше- 
му, а при одинаковых суммах — обычное лексикографическое упорядочение слева направо от большего к меньшему; 


а. 
А ИХ 


м] — КОЛЬЦо Полиномов степени не выше отт переменных Над полем В : 


Степень 4е2(/) полинома Г: определяется как максимальная степень составляющих его ненулевых моно- 


мов, а степень ненулевого монома ф=х\...х” =ах" задается равенством 4е5(ф) =р(у). 


“би 


— Г 2 и 
Вектор УЕР№, где К=|+т+С»..1, составленный из коэффициентов информационного полинома 


ие ® 
У(х) = Е нь Хх ]. называется информационным вектором. При этом предполагается, что для нумерации эле- 


ментов информационного вектора, как и кодового, используется упорядочение (1). Кодирование информационного 


РИД 


вектора осуществляется путем вычисления значений соответствующего информационного полинома в точках про- 
странства у ‚ упорядоченного в соответствии с (1). Способы нахождения числа гарантированно исправляемых оши- 
бок в общем виде см. в [4], [8]. Далее в работе будем использовать уже вычисленные значения. 

Конструкция ДМР для кодов МЗ (2,т). На вход алгоритма подаются параметр т кода АМз(2,т), связанные с т зна- 
чения длины и и размерности К кода ЮМз(2,т), а также полученный из канала зашумленный кодовый вектор 
У = (У, .....7, ЕЕ. (= С”), элементы которого занумерованы в соответствии с (1). На выходе алгоритма формируется 


восстановленный информационный вектор у ; 


, р И Я Е —1 —1 
Шаг 1. Построим набор векторов из С”: У.) = (а, а а , 
Г. у 


Ти 2 +0, 0 2 ^ ЕРу ‚+0 





где у. ' — число, сопряженное У, ‚ векторы 7 е Ру’ упорядочены по (1). Далее везде, где нумерация элементов векто- 


р т х 
ров или наооров осуществляется векторными переменными из Е ‚ По умолчанию оудем использовать упорядочение 


(1). 


ПТаг 2. Рассмотрим все значения УЕ К”, у#0 и для фиксированного 7 введем обозначение 


и 
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В=(В, =0,В, „..,В.) И А=(А, =0,4, ,...А,) Б, где В, — вектор 


В = (Р......В„,)е Ру’, на котором достигается А. — минимальное значение функционала 


АБВ =У” В, 69" (с К), 


0. 





где (Ва, )( №) — скалярное произведение, В, ЕЁ.. 
Шаг 3. Построим (Их) -массив @, строки которого инициализируем значениями из набора В: 


©(а,)=В. =(0,. .....0,=.)е Е”. Далее 7-Й столбец полученного массива будем обозначать (-) у = 1..... т. Для каж- 


я 


дого 09, ЕР; и всех В, =Ру” таких, что В, = ©, ВЫЧИСЛИМ 


@(а,, ) = Ма (@(а, +В,)-©(В,)), нра, 


где функция Ма возвращает элемент, встречающийся наибольшее число раз. 
Шаг 4. Для каждого 1 =]...., 2 найдем Я ‚ как минимум функционала 


я 


Г, (ф) = Е Аь, 





(ЕК), 


7 





заданного на множестве всех линейных однородных полиномов вида ф(х)= р ф м ф. ЕЁ. Полином ф, на кото- 


4=1 


— ие т га) 

5ом достигается минимум. обозначим ®7’(х)= ох. 
Е = 9 Ч 
х = 


Вычислим =>. Ча Х) ‚где хе В", 4, ЕВ, 4,] Е [1..,т], а, =а, и 


(а) : 

р @, ‚если |} и. Я , 
7 |007 еслиа. за 
ат та / 


Шаг 5. Среди множества векторов 6 = (6,.....С, Е В” и значений С, ЕД найдем те, которые минимизируют 


$7,5) => 





у, —_ ео Со +55, )-+ ча, )} К К) 


где ие 0. ) (= Ё, ) — скалярное произведение. 


“^ `„ -—- 7 
Из найденных значений составим полином ф(х) = ут „х, < ‚ где коэффициенты су © Рз, А=О,.., т со- 


ответствуют найденным значениям С И 5 р 


Результат декодирования строим в виде полинома: /(х)=\и/(х)-+Ф(х), который определяет искомый инфор- 


«- 


мационный вектор у ь 


Экспериментальное исследование. В работе проведено моделирование передачи закодированной информации по 
троичному каналу связи как с жесткими, так и с мягкими решениями приемника. Канал связи моделировался с ис- 
пользованием информационной системы «Канал» [9], для которой были созданы специальные библиотеки. При про- 
ведении экспериментов параметры модели задавались следующими входными данными: 

— значение т, определяющее параметры помехоустойчивого кода АМз(2,т); 

— число ошибок +, поражающих кодовое слово; 

— тип приемника. 

Если использовался приемник с мягкими решениями, то применялись дополнительные настройки, указываю- 
щие на вид используемых базовых искажений элементов кодовых слов, а также параметр приемника 5, задающий 
допустимую область значений =, , при использовании искажений по амплитуде. Проведено 10“ испытаний для каждо- 
го набора параметров модели. 

Отметим, что при моделировании потоков ошибок часто используется понятие пораженных символов, т. е. 
попавитих под воздействие искажений, но по случайности не изменивших своего значения [3], [10]. При проведении 
экспериментов такая ситуация намеренно не моделировалась, т. е. в результате наложения # ошибок на кодовое слово 
это слово отличалось от исходного ровно в Е позициях. 
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В табл. | описаны основные параметры кодов Ю/Мз(2,т), результаты экспериментального исследования кото- 
рых представлены ниже. В этой таблице использованы следующие обозначения: т — параметр кода ЮМз(2,т); п, Кит 





— длина, размерность кода и число гарантированно исправляемых ошибок (определяемое по минимальному расстоя- 
нию кода). Отношение и задает максимальное значение вероятности ошибки в канале связи, при котором код «гаран- 
тирует» выдачу верного результата. Отметим, что при ошибках, вероятность которых превышает и, классический 
детерминированный декодер всегда выдает ошибочные результаты. Параметр Ии показывает избыточность кода. При 
необходимости далее будем использовать традиционную краткую запись параметров кода в форме тройки [и, А, 4], где 
4 — минимальное кодовое расстояние. 

Таблица | 

Основные параметры кодов АМз(2,т), т=2, 3, 4, 5 


ыд 
О МС ПО ООС ООО 





Проведенная серия экспериментов по определению корректирующей способности нового ДМР-кода АМз(2,т) 
в случае применения его в канале с жесткими решениями приемника показала значительное повышение числа ис- 
правляемых ошибок по сравнению с детерминированным декодером. Изменение значения максимальной вероятности 
исправляемых оптибок новым декодером по сравнению с детерминированным декодером представлено в табл. 2. В 
верхней строке таблицы указаны максимальные вероятности ошибок, исправление которых гарантируется кодом (па- 
раметр #/и, см. табл. 1). Вторая строка таблицы содержит максимальные вероятности ошибок, при которых новый мяг- 
кий декодер выдал верный результат во всех проведенных экспериментах (т. е. в 100 % случаев). Из приведенных ре- 
зультатов видно, что декодер значительно улучиил результат, гарантированный кодом, во всех случаях, кроме кода 
КМ3з(2,2). Так для КМз(2,3) корректирующая способность увеличилась на 34 %, для АМз(2,4) — на 54 %, а для КМ3з(2,5) 
— на 157 % (код гарантирует исправление 40 ошибок на кодовое слово, а декодер во всех проведенных экспериментах 
исправил все ошибки до 103 включительно на кодовое слово). Гретья и четвертая строки таблицы отличаются от вто- 
рой вероятностью выдачи декодером верного результата (третья строка — верный результат в 95 % случаев, четвертая 
— в 90 % случаев). 
Таблица 2 
Корректирующая способность нового декодера кода АМ3з(2, 171 ) 


в случае применения его в канале с жесткими решениями приемника 


Максимальное значение вероятности исправляемых Параметр кода Е т) 


кодом АМз(2,т) т 11 0,160 0,165 


декодером с вероятностью 1 0,111 0,148 0,247 0,424 
декодером с вероятностью 0,95 0,111 0,185 0,321 0,428 
декодером с вероятностью 0,9 0,111 0,185 0,333 0,436 





Для определения чувствительности нового ДМР-кода АМз(2,т) к различным базовым типам искажений и их 
комбинациям проведена серия экспериментов с различными настройками параметров модели канала связи. На рис. 4 
представлен график зависимости вероятности верного декодирования от вероятности ошибки в канале связи для кода 
ЮМ3з(2,4). 

Рассмотрим график подробнее. Кривая с подписью «Дискретные ошибки» отражает результаты эксперимен- 
тов, проведенных в канале связи с жесткими решениями приемника. Напомним, что в этом случае ошибка переводит 
сигнал в один из двух других разрешенных сигналов из поля Сз. Кривая с подписью «Ошибки по фазе» отражает ре- 
зультаты экспериментов, проведенных в канале связи, где происходят только ошибки типа сдвига сигнала по фазе. 
Кривые с подписью «По фазе и амплитуде» отражают результаты экспериментов, проведенных в канале связи, где 


происходят различные комбинации ошибок по фазе и по амплитуде. в скобках указан Параметр приемника 5, опреде- 
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рии = 


ляющий ширину допустимой области =.. Легко видеть, что чем меньше значение &, тем область =, больше. Экспе- 
рименты показали, что чем итире область допустимых значений (Т. е. чем значения & меньше), тем хуже корректиру- 
ющая способность декодера. Однако для каждого кода находится такое значение &, при увеличении которого значи- 
тельного изменения корректирующей способности практически не происходит. Гак для кода АМз(2,4) результаты экс- 
периментов при значениях &Е [0.5, |] практически неразличимы. На графике не представлена зависимость корректи- 
рующей способности декодера в случае, когда происходят только опгибки по амплитуде. Эксперименты показали, что 
декодер не чувствителен к таким ошибкам, если не происходит сдвига сигнала по фазе. Так, для числа ошибок типа по 
амплитуде от | до длины И кода декодер всегда выдает верный результат. Отметим, что взаимное положение кривых, 


представленных на графике, повторяется для различных кодов АМз(2,т), т>2. 


=== [10 фазе и амплитуде (0.9) 
=ф== [По фазе и амплитуде (0,5) 
А = == []о фазе и амплитуде (0.1) 


-=------=---- р нннн----- 
, 


Еее: 


== Дискретные ошибки 


Вероятность верного декодирования 





0,1 0,3 0,5 0,7 
Вероятность ошибки в канале 


Рис. 4. Результаты исследования [81, 15, 27]|-кода К Мз(2.,4) 


В работе проведена серия экспериментов по поиску зависимости корректирующей способности декодера от 
расположения ошибок внутри кодового слова. Результаты экспериментов показали, что декодер не чувствителен к 
местоположению ошибок. 

Из табл. [| видно, что [243, 21, 81 |-код ЕМ3(2.5) обладает достаточно высокой избыточностью и/к = 11,6. Воз- 
никает очевидный вопрос: целесообразно ли использовать этот код с его трудоемким алгоритмом декодирования или 
имеет смысл использовать коды с болышой избыточностью и простыми декодерами? Рассмотрим в качестве кода, 
имеющего хорошую корректирующую способность и простой алгоритм мажоритарного декодирования, код повторе- 
ния |1], [11] с близким значением избыточности и сравним корректирующие способности кода ЮМз(2.5) и [11,1,11]- 
кода многократного повторения над алфавитом Ёз, который гарантированно исправляет 5 ошибок на | кодовое слово 
и обладает избыточностью, равной 11. Для того, чтобы сравнять длины кодовых слов обоих кодов, объединим 22 ко- 
довых слова кода повторения — например, с помощью техники прямой суммы кодов [1], [11]. Гаким образом, длина 
составного кодового слова — 242 символа, и код может гарантировано исправить [10 символов в этом составном сло- 
ве, но только при жестком ограничении: ошибочные символы должны быть распределены внутри кодового слова рав- 
номерно (не более 5 ошибок на 1]-символьный отрезок кодового слова). Код Рида — Маллера может гарантированно 
исправить только 40 ошибок внутри кодового слова длиной 234 символа, однако результат его работы не зависит от 
местоположения ошибок. Если использовать новый мягкий декодер кодов Рида — Маллера, то в рассматриваемом 
случае можно будет исправить 103 ошибки, также независимо от места их положения. Более того, если применять 
новый декодер в полунепрерывном канале связи, то декодер сможет исправить не менее 150 ошибок, вне зависимости 
от места их положения внутри кодового слова. Декодер кода многократного повторения не предназначен для исполь- 
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зования в полунепрерывном канале связи. Очевидно, что код многократного повторения, который обычно рассматри- 
вают как код, обладающий большой корректирующей способностью, уступает по этому параметру коду ЮМ3з(2,5). хо- 
тя, очевидно, выигрывает по скорости работы алгоритмов кодирования и декодирования. Для кодов АМз(2,т) с дру- 
гими значениями параметра т аналогичные рассуждения также показывают выигрыш в корректирующей способности 
по сравнению с кодом многократного повторения. При этом с ростом т данный выигрыш растет, а с уменьшение зна- 
чения 721 — падает. 

Заключение. Результаты проведенных экспериментов показали, что исследуемый декодер троичного кода Рида — 
Маллера второго порядка обладает значительной корректирующей способностью по сравнению с классическими де- 
терминированными декодерами — например, с декодированием по минимуму расстояния Хэмминга. Новый декодер 
мягких решений может быть применен для обеспечения помехоустойчивости в каналах связи низкого качества, ис- 
пользуемых, однако, для передачи ценных сообщений. Так, на графике 4 видно, что даже при вероятности ошибки в 
канале связи, равной 0,5, т. е. значению, при котором говорят о разрыве линии связи, декодер выдает верные результа- 
ты с высокой достоверностью. Кроме традиционного применения декодера для обеспечения помехоустойчивости он 
может быть использован в таких практических приложениях, как восстановление информации, полученной по побоч- 
ным (отводным) каналам связи [12], при построении криптосистем типа Мак-Элиса и Ниддерайтера [11]. 
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